Draft version July 23, 2010 

Preprint typeset using I^'T^]X style emulatcapj v. 2/16/10 



A CORRELATION BETWEEN STELLAR ACTIVITY AND HOT JUPITER EMISSION SPECTRA^ 
Heather A. Knutson^-^, Andrew W. Howard^-^, & Howard Isaacson^ 

Draft version July 23, 2010 



o 

(N 



(N 
(N 

o 



(N 
> 

O 
l> 
(N 

^' 

O 
O 



ABSTRACT 

We present evidence for a correlation between the observed properties of hot Jupiter emission 
spectra and the activity levels of the host stars measured using Ca H H & K emission lines. We find 
that planets with dayside emission spectra that are well-described by standard ID atmosphere models 
with water in absorption (HD 189733, TrES-1, TrES-3, WASP-4) orbit chromospherically active stars, 
while planets with emission spectra that are consistent with the presence of a strong high-altitude 
temperature inversion and water in emission orbit quieter stars. We estimate that active G and K stars 
have Lyman a fluxes that are typically a factor of 4 — 7 times higher than quiet stars with analogous 
spectral types, and propose that the increased UV flux received by planets orbiting active stars destroys 
the compounds responsible for the formation of the observed temperature inversions. In this paper we 
also derive a model-independent method for differentiating between these two atmosphere types using 
the secondary eclipse depths measured in the 3.6 and 4.5 /im bands on the Spitzer Space Telescope, 
and argue that the observed correlation is independent of the invcrted/non-inverted paradigm for 
classifying hot Jupiter atmospheres. 

Subject headings: binaries: eclipsing — stars: activity — planetary systems — techniques: spectro- 
scopic 



1. INTRODUCTION 

The close-in, gas giant planets known as "hot Jupiters" 
occupy a unique regime of parameter space, with tem- 
peratures between 1000 — 2500 K and atmospheric com- 
positions that are likely to be quite different than that 
of the cooler solar system gas giants. Because the 
time scale for tidal synchronization is short compared 
to the ages of these systems, hot Jupiters in circular or- 
bits are expected to be tidally locked, leaving one side 
of the planet in permanent darkness while the other 
is constantly illuminated by intense radiation from the 
star. This further complicates the atmospheric chem- 
istry, as the cooler night side may act as either a 
source or a sink for different materials, while vigor- 
ous day-night circulation may produce increased verti- 
cal mixing on the dayside. These planets also experi- 
ence much higher UV fluxes than the solar system gas 
giants, and it is likely that photochemistry alters the at- 
mosp heric chem istry down to pressures of 10 mbars or 
more (iLiang et al. 2003, 200l iZahnle et all [2001 120T0I: 
iLine. Liang, fc Yuna i2010i) . Infrared light from these 
planets emerg es fro m pressures of 10-100 mbar (e.g. 
iBurrows et al.l l2008': Sh owman et ani2009l ). and as a re- 
sult their dayside emission spectra could exhibit strong 
features from photochemically produced molecules. In- 
deed, the detection of an excess o f CO2 in the atmosp here 
of the hot Jupiter HD 189733b (jSwain et al.ll2009aD in- 
dicates that non-equilibrium chemistry may play an im- 
portant role in these atmospheres. 

^ Based on observations obtained at the W. M. Keck Observa- 
tory, which is operated jointly by the University of Cahfornia 
and the California Institute of Technology. Keck time has been 
granted by both NASA and the University of California. 

^ Department of Astronomy, University of California, Berke- 
ley, CA 94720-3411, USA 

^ Miller Fellow; hknutson@berkeley.edu 

* Townes Fellow, Space Sciences Laboratory, University of Cal- 
ifornia, Berkeley CA 94720-7450 USA 



We can characterize the dayside emission spectra of 
these planets by measuring the wavelength-dependent 
decrease in light as the planet passes behind the star 
in an event known as a secondary eclipse. Secondary 
eclipse observations have now been carried out for more 
than sixteen extrasolar planets, and there is currently a 
growing body of evidence for the existence of two d i stinct 
classes of hot Jupiter atmospheres (see iDemind I2009L 
and references therein). According to the current clas- 
sification scheme one cla ss of planets (e.g. HD 189733; 
iCharbonneau et aI1l2008D have dayside emission spectra 
that are well-described by atmosphere models with wa- 
ter and CO in absorption . A second class of planets (e.g. 
HD 209458 : [Knutson et a l. 2008) are better described by 
models with a strong temperature inversion between 0.1- 
0.01 bars and these same features in emission. 

We do not know the nature of the absorber re- 
quired to create these inversions; it was initially sug- 
gested that gas-phase TiO might provide the nec essar; 
opacitvdHu bcnv et al.' '2003t IBurrows et all 120071 1200 
iFortnev et al. 2006 , 20081 ) . but this theory cannot ex- 
plain the full range of current observations. TrES-3 is 
hot enough for g a.s-phase TIO but d oes not appear to 
have an inversion (jFressin et al.l[20Toh . whereas XO-1 re- 
ceives a much lower incident flux comparable to t hat of 
HD 1 89733b and yet has an inversion (Machale k et al.l 
|2008[ ). It is also unclear whether or not TiO can be reli- 
ably maintained in the upper atmosphere at all, given its 
high molecular weight and the likely presence of day- and 
night-side cold traps. Overcoming these effects would 
require vigorous mixing, but this niay not be expecte d 
in a stably stratif ied atmosphere (iSpiegel et akl 12009) . 
More recently, iZahnle et al.l (|2009[ ) presented detailed 
non-equilibrium atmospheric chemistry models suggest- 
ing that heating from sulfur compounds in the upper at- 
mospheres of hot Jupiters could explain these inversions; 
in this model, photochemistry serves as a net sink for the 
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Figure 1. Ca II H & K lines for transiting planet host stars ob- 
served with Keck/HIRES, plotted in order of increasing Shk- This 
plot shows half of our sample, with the second half plotted in Fig. 
[2] Active stars have significant emission in the line cores (high val- 
ues of ^^'^ li'' the top of the plot; a majority of the stars in 
our sample appear to have no detectible emission. Star names are 
given on the left, with systems hosting planets with temperature 
inversions shown in red and those without temperature inversions 
shown in blue, while names in black denote systems with insuffi- 
cient data. It is possible that cooler (< 1000 K) and/or higher 
density planets, such as the Neptune-mass planets GJ 436 and 
HAT-P-11, may not fit into the simple classification scheme used 
to describe hot Jupiters atmospfieres (indeed this appears to be 
the case for GJ 436; Steven son et al.ll20ig) . but we include HIRES 
observations in this plot where available for completeness. 



sulfur compounds of interest. 

In order to quantify the effects of photochemistry 
in these atmospheres more precisely, we must first ob- 
tain an estimate of the incident UV flux. Direct 
measurements of the stellar UV fluxes are only avail- 
able for a handful of systems, in cluding HD 189733b 
(iLccavclier des Etangs et all l2010f) and HD 20945813 
■■Vidal-Madiar et al. 2003'. '2004: 'Ehrenr eich eral| [205l 



JVi 



France et al. 2010; Linsky et al. 2010), while arc hival 
searches for X-ray detections ([Kashvap et al.l 120081 ) and 
more recent ob servations of planet-hosti ng stars with 
XMM-Newton (jPoppenhaeger et al.ll2010l ) provide only 
upper limits in most cases. In this paper we prese nt mea- 
surements of Ca II H & K line strengths (,Noves et al.l 
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Figure 2. Ga II H & K lines for transiting planet host stars 
observed with Keck/HIRES, plotted in order of increasing 5jjk- 
This plot is a continuation of Fig. [T] and includes stars with lower 
values of ^jj^ • 



Il984( ). which act as an indicator of stellar activity lev- 
els, based on Keck HIRES spectra for fifty transiting 
planet host stars. These observations allow us to iden- 
tify the most active planet-hosting stars, which should 

have correspondingly enhanced UV and X-ray fluxes 
^„ ... „ . „ „ .. . 



Schrijver. Dobson. fc Radic3ll992t [Liv ingsto n et al.l 



12007 ). and to search for evidence that increased activity 
in the host star can explain the observed dichotomy in 
hot Jupiter emission spectra. We describe our observa- 
tions in |J2] and survey the available literature on clas- 
siflcations of hot Jupiter emission spectra in ^J3l In ^J4] 
wc develop a model-independent metric for distinguish- 
ing between the two observed types of planetary atmo- 
spheres, and argue that planets without inversions are 
consistently found around the most active stars, while 
planets with strong inversions orbit quiet stars. We then 
estimate the effect that increased activity has on the mea- 
sured Lyman a fluxes from the stars in our sample in ij5l 
and discuss the consequences this has for our understand- 
ing of hot Jupiter atmospheres in fJSl 

2. KEGK/HIRES Ga II H & K OBSERVATIONS 



A Correlation Between Stellar Activity and Hot Jupiter Emission Spectra 
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Figure 3. Ca II H & K lines for a sample of seven F stars observed 
with Keck/HIRES as part of ongoing CPS programs, plotted in 
order of increasing 5jjj^ (values given on left). These stars were 
selected to have effective temperatures between 6000 — 6500 K, 
vsini less than 10 km s~^, and ^jj^ values in the range 0.18 — 0.24, 
comparable to the values measured for XO-3 (S'jjj^=0.239) and HD 
147506 (Sjjx^O-lSl)- The original spectra (thin grey lines) for 
these stars were then artificially broadened to a v sin i of 20 km 
s~^ (thick black lines), comparable to that of XO-3 (second from 
top) and HD 147506 (second from bottom) which arc plotted for 
comparison. The effective temperatures for the stars shown in this 
plot are, from bottom to top, 6170, 6080, 6260, 6090, 6130, 6090, 
6190, 6430, and 6090 K, respectively 

We observed the transiting planet host st ars listed in 
Table[T]with the HIRES echelle spectrometer ()Vogt et alj 
I1994D on the 10-m Keck I telescope (see Fig. [T]and[2l). 
We used the standard HIRES setup e mployed by the Cal- 
ifornia Pla net Search (CPS) group (jWright et all 120041: 
[Howard "et* al. 2009). Most observations were made with 
a slit of width 0'.'86 for 500 s yielding a signal to noise of 
75-200 per pixel (depending on the host star brightness) . 
Many of the observations were initially made for radial 
velocity monitoring, and additional single-exposure ob- 
servations were obtained where necessary to expand this 
coverage to a larger sample of transiting planet host 
stars. For the radial velocity observations, an iodine 
cell was mounted directly in front of the spectrometer 
entrance slit to provide a wavele ngth scale and calibra- 
tion of the instrum ental profile ()Marcv fc Butleii Il992l : 
iValenti et al.lfT995l ). The iodine absorption lines (5000- 
6200 A) do not affect the measurement of Ca II H & 
K indices (3968 A and 3933 A). Stellar parameters were 
determined using a LTE spectral synthesis routine (SME 
or Spectroscopy Made Easy; Fischer & Valcnti 2005) on 
observations taken without the iodine cell where avail- 
able. Systems with only a single spectrum were observed 
in this mode. 

We calibrate the stellar activity indices to the Mt. Wil- 
son S'-value scale, defined as the ratio of the sum of 
the flux in the cores of the Ca II H & K lines to the 
sum of two continuum bands, one redward a nd one blue- 
ward of the H and K lines (|Wilson I IT968I ) . The cali- 
bration from flux to Suk was accomplished using stars 
that were observed i n both the Mt. Wilson H-K project 
(jPuncan et al.|[r991 ') and the California Planet Search. 
Shk values for a particular star on the CPS scale are typ- 



ically within 10% of the Mt. Wilson scale (jlsaacson et al.l 
120101 ). Before extracting the H hne at 3968.5 A, the K 
line at 3933.7 A and continuum flux levels, the spectra 
are shifted to rest wavelengths through a cross correla- 
tion of the NSO solar atlas. The flux in the H and K lines 
cores is measured in two 1 A FWHM weighted triangles 
centered on the H and K lines. The continuum sections 
are each 20 A wide and reside sufficiently far away from 
the H and K line centers so as to avoid the broad wings 
of either line. After calculating the S- value for each star, 
we convert to \o^(R '^^) using the B — V color of the star 
(|Noves et al.l 1198^ : this effectively removes any depen- 
dencies on the bolometric flux of the star and allows us 
to accurately compare the chromospheric emission from 
different spectral types. In cases where we have multiple 
spectra available spanning several epochs, we take the 
median value of log(i?^j^) over all of the available obser- 
vations. We list the resulting Shk and log(i?^j^) values 
in Table [TJ for a more detail ed description of our o bser- 
vations and methodology see llsaacson et alJ (|2010( ) . 

2.1. Shk o,nd log(R'jj^) for F Stars 

Although we calculate S'hk and log(i?yj^) values for the 
F stars (Tofr=6000 — 6500 K) in our sample, we note that 
the lo g(i?lrT^) scaling relations described in iNoves et al.l 
(|1984| ) are not well calibrated for B — V< 0.5. Because 
F stars have higher continuum fluxes in the region of 
the Ca II H & K lines, it can be difficult to detect 
small amounts of emission in the line cores even when 
it is present. Our values for S'hk and log(i?HK) -^O" 
3 (Toff=6430 K, t;sini=18.5 km s'^) and HD 147506 
(Toff=6260 K, wsini=20.7 km s~'^) suggest that both 
stars are active, but these values should be regarded with 
some suspicion as these stars are on the edge of the cal- 
ibrated range in B — V. 

As an additional test, we carried out a visual inspec- 
tion of the spectra for XO-3 and HD 147506 and could 
find no evidence for emission in the Ca II H & K line 
cores. It is possible that the rotational broadening could 
be obscuring the presence of weak lines; we test this the- 
ory by searching the database of all CPS spectra and 
selecting a sample of F stars with effective temperatures 
and S'hk values similar to those of XO-3 and HD 147506 
(6000 < Teff < 6500 K, 0.18 < 5hk < 0.24), but 
with significantly lower (< 7 km s^^) values of vsini. 
We then artificially broaden these spectra to a w sin i of 
20 km s"\ comparable to that of XO-3 and HD 147506, 
using the standard IDL routine lsf_rotate.pro, available 
as par t of the astrolib distribution and described in lGravl 
(|2005f) . A characteristic double-peaked emission feature 
is clearly detectible in the line cores of the unbroadened 
spectra (see Fig. |3|), and although this signal is obscured 
somewhat by the broadening, there is still a noticeable 
flattening or even an increase in flux near the centers of 
the line cores in the high wsini versions of the spectra. 
In contrast to these systems, XO-3 and HD 147506 (also 
plotted in Fig. [3] for comparison) both have v-shaped 
spectra more characteristic of quiet stars. We conclude 
that both stars are likely to be chromospherically quiet, 
with the additional caveat that the standard Shk and 
log(i?^j^) indices used to measure this emission are not 
well-cahbrated and are probably not appropriate to use 
for stars of this spectral type. 
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Table 1 
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& K Line Streng 






Star 


T,s° (K) 


V sini'' 


B-V 






HD 17156 


5980 


4.3 


0.64 


0.155 


-5.022 


HD 147506 


6260 


20.7 


0.41 


0.191= 


-4.780 


HD 149026 


6150^= 


6.2 


0.61 


0.152 


-5.030 


HD 179949 


6170 


7.0 


0.50 


0.232 


-4.622 


HD 189733 


5090 


2.7 


0.93 


0.508 


-4.501 


HD 209458 


6070 


4.0 


0.59 


0.160 


-4.970 


xlU oUDUD 


OOiU 


2.4 


77 


U. 104 


-O.UDi 


GJ 436 


359QC 


< 3 


1.52 


0.620 


-5.298 


V And 


6160^= 


9.7 


0.54 


0.156 


-4.982 


TrES-1 


5300 


0.5 


0.78 


0.244 


-4.738 


TrES-2 


5840 


1.6 


0.62 


0.165 


-4.949 


TrES-3 


5530 


3.1 


0.71 


0.305 


-4.549 


TrES-4 


6200== 


8.5 


0.52 


0.141 


-5.104 


XO-1 


5750= 


1.1 


0.69 


0.168 


-4.958 


XO-2 


5340= 


1.4 


0.82 


0.173 


-4.988 


XO-3 


6430= 


18.5 


0.46<^ 


0.239= 


-4.595 


XO-4 


6400= 


8.8 


0.47"^ 


0.124 


-5.292 


HAT-P-1 


5980 


2.2 


0.58<^ 


0.158 


-4.984 


HAT-P-3 


5170 


0.5 


0.87<* 


0.206 


-4.904 


HAT-P-4 


5990 


5.1 


0.58<^ 


0.145 


-5.082 


HAT-P-5 


5960= 


2.6 


0.59<i 


0.148 


-5.061 


HAT-P-6 


6410 


8.5 


0.41 


0.187 


-4.799 


HAT-P-7 


6350 


3.8 


0.44 


0.150 


-5.018 


HAT-P-8 


6130 


11.5 


0.54^ 


0.156 


-4.985 


HAT-P-9 


6350= 


11.9 


0.48'* 


0.141 


-5.092 


HAT-P-IO'' 


4990 


0.5 


1.01 


0.322 


-4.823 


HAT-P-11 


4820 


0.5 


1.02 


0.580 


-4.567 


HAT-P-12 


4650 


0.5 


1.13<* 


0.253 


-5.104 


HAT-P-13 


5710 


2.7 


0.73 


0.141 


-5.138 


HAT-P-14 


6530 


8.5 


0.42 


0.175 


-4.855 


HAT-P-15 


5570= 


2.0 


0.71^ 


0.166 


-4.977 


HAT-P-16 


6160= 


3.5 


0.53<* 


0.175 


-4.863 


WASP-1 


6170 


0.2 


0.53d 


0.140 


-5.114 


WASP-2 


5230 


1.3 


0.84<* 


0.159 


-5.054 


WASP-3 


6170 


14.2 


0.52 


0.173 


-4.872 


WASP-4 


5500= 


2.0 


0.74'* 


0.194 


-4.865 


WASP-12 


6300 


3.0 


0.50'* 


0.113 


-5.500 


WASP-13 


5910 


5.0 


0.60'* 


0.127 


-5.263 


WASP- 14 


6270 


3.9 


0.45 


0.163 


-4.923 


WASP-17 


6380 


iU. / 


n A fid 
U.4o 


0.121 


-5.331 




6250 


11.7 


0.49 


0.116 


-5.430 




ooyu 


5.3 


0.70'* 


u. zoz 


-4.DDU 


CoRoT-1 


5950= 


5.2 


0.59'* 


0.124 


-5.312 


CoRoT-2 


5630= 


11.5 


0.69'* 


0.435 


-4.331 


CoRoT-7 


5330 


2.5 


0.80'* 


0.225 


-4.802 


Kepler-4 


5860= 


2.2 


0.62'* 


0.168 


-4.936 


Kepler- 5 


6300= 


4.0 


0.50'* 


0.148 


-5.037 


Kepler-6 


5650= 


3.0 


0.68'* 


0.160 


-5.005 


Kepler-7 


5930= 


4.2 


0.59'* 


0.155 


-5.099 


Kepler-8 


6210= 


10.5 


0.52'* 


0.153 


-5.050 



Also known as WASP-1 1. 
^ Values for Tefi (K) and vsini (km s^*) were determined from 
a SME analysis of th e Keck/HIRES spectra as described in 
IFischer &: ValentH l|2005l ') where available, otherwise values from the 
literature were used instead. 

= We revert to the published values for stellar effective temperature 
and tisini in these cases. 

=* In some cases measured B~V values were unavailable or were 
inconsistent with the star's spectral type. In these cases we opted 
to set the B — V col or for each star equal to the value given in Table 
B.l in I Gray! (120051 ) for a given effective temperature. We do not 
account for variations in B — V due to metallicity or age, as these 
quantities have a negligible effect on our final log(Hj^j^) values for 
the range of stellar spectral types present in our sample. 
= These hot, rapidly rotating F stars have Sjjk ^.nd log(_R^j^) val- 
ues that are suggestive of activity but they are near the edge of 
the range hi B — V for which these indices are calibrated, and their 
spectra show no detectible emission in the Ca II H & K line cores. 
We argue in ^2.1l that these stars are most likely chromospherically 
quiet. 

* Shk 3.nd log(_R^j^) values are only calibrated for stars with 
B-V between 0.5. - 1.4 (4200 - 6200 K). We give values for all of 
the stars in our sample, but values for spectral types outside this 
range should be treated with some skepticism. 



3. OBSERVATIONS OF HOT JUPITER EMISSION SPECTRA 

There are currently four planets classified in the lit- 
erature as having "non-inverted" atmospheres, includ- 
ing HP 189 7 33b flPeming et al 



20071 120081 : iCharbonneau et al 



2006t iGril lmair et"aT 



12008: Swain et al 



2009al). TrES-1 (ICharbonneau et al.i i2005ll. T rES-3 
(jFressin et all 120101 ) . and WASP-4 (|Beerer et all [20101. 
These planets have emission spectra that are consis- 
tent with standard ID cloud-free atmosphere mod- 
els where temperature de creases with increasing height 
in the atmosphere (e.g . 'Barman' '2 0081: iBurrows et alj 
[2008l:lFortnev et al.ll2Q08 : Madhusudh an fc Seaged2009ir 
or display at most a weak inversion (e.g. WASP- 
4; iBeerer et al.l l2010l ). Models used to describe these 
planets exhibit CO and H2O features in absorption, 
which provides a good match to the shape of their ob- 
served emission spectra in the mid-IR (3.6-24 /im). Ob- 
servations of TrES-1, TrES-3, and WASP-4 are lim- 
ited to a subset of the Spitzer IRAC bandpasses (3.6, 
4.5, 5.8, 8.0 /im), which have relatively wide (0.75 — 
2.9 /im FWHM) transmission functions and are there- 
fore somewhat ambiguous in their interpretation, but 
HD 189733b has been observed at higher wavelength 
resolutions with bo t h IR S on Spitzer (5 — 14 /im; 
IGrillmair et"aLl [20071 [200l and N ICMOS on the Hub- 
ble Space Telescope (1.5 — 2.5 /im: iSwain et a l. 2009a). 
Published models for HD 189733b in which the at- 
mospheric chemistry and pressure-temperature profiles 
are allowed to vary in the fit consistently find that 
this planet is best-described by a pressure-temperature 
profile that decreases with height in the atmosphere 
([Madhusudhan fc S eage? 2009; Swain et al."2009a'). 

In contrast, the emiss ion spectra o f planets such as 
HD 209458b (Domin g etaU [2005t [Richardson et"all 
[2007t iKnutson et al.. ,2008t ISwain et al.l I2009bl). TrES- 
2 ([O ^Donovan et al.l I2010D. TrES - 4 (IKnutson et alJ 
I2009al). HAT-P-1 (Todorov et all I2010D. HAT-P-7 
(iChristiansen et al.l .2010), W ASP-lbfWheatlev etall 
I2010D. WASP- 18 b (pfvmever et al.l I2010D. XO^ 
(iMachalek et al.l l200 8[), X0-2b ([Machalek eFall 
I2009D. X0-3"bl Machaiek et al.l [20M). and CoRoT-lb 
([Deming et al.l [2010, 1 are best-described by models with 
a stronger temperature inversion betwee n 0.1-0.01 bar 
and water and CO bands in emission (e.g. Hubenv et aL[ 
200l iFortnev et all [2001 [200l [filirrows et al.. .200 



2008; iMadh usudha n fc Seaged 120091: iSpiegel fc Burrowsl 
2010[) . As before, our information on the majority 
of these systems is limited to observations in the 
Spitzer IRAC bandpasses, although HD 209458b 
has also been obser ved with IR S (Ric hardson et alJ 
[20071 ) and NICMOS ([Swain et al.l [2009b1 and shortcr- 
wavelengt h data is available for a few add i tional planets 
jBorucki c iTall [20091: lAlonso et all [20091: [Snellen et af 
2009t iChristiansen et all 120101 : iLopez-Morales et al 
20101) . Model fits to HD 209458b's emission spectrum in 



which the chemistry and pressure-temperature profiles 
are allowed to vary find that this planet is best-described 
by a model with a ternperature inversion at al titude 
([Madhusudhan fc Seage3l2009l: ISwain et al.ll2009bO . con- 
sistent with the conclusions of m ore highly-constrained 
models (e.g. Burrows et al."2007'). 

Ahhough Wheatlcy ct al. (201(| classify WASP-2b as 
a non- inverted planet, this interpretation relies on a spe- 
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Figure 4. Stellar effective temperature vs. Ca II H & K activ- 
ity index log{R'^^), which tracks the amount of emission in the 
Ca II H & K line core. Active stars have stronger emission fea- 
tures and correspondingly large values of log(i?yj^); because the 
calibration for log(i?^j^) is uncertain for stars with effective tem- 
peratures higher than 6200 K IINoves et al.lll98^ , we have marked 
this region with diagonal grey lines on this plot. Planets classified 
in the literature as having temperature inversions are shown as 
red circles, while planets in the literature that are well-described 
by non-inverted atmosphere models are shown as blue stars (see 
U for a discussion of WASP-2). CoRoT-2 (grey square) has an 
unusual spectrum that is not well-described by either model, and 
XO-3 (red downward arrow) has temperature inversion and orbits 
a star with a relatively high log(ij^j^) value but no visible emission 
in the Ca II H &; K line cores by inspection (see J2TT}. Although 
cooler stars tend on average to be more active, there is still a range 
of activity possible within a given spectral type: XO-2, WASP-2, 
and TrES-1 differ by less than 100 K, and yet TrES-1 is moderately 
active while XO-2 and WASP-2 appear to be quiet. For reference 
we give the spectral types for a range of effective temperatures at 
the top of the plot. 

cific model for creating the inversion (a solar metallicity 
atmosphere and gas phase TiO in local thermal equili- 
brum) that does not always provide a good match for 
the observed spectr a of other planets (e.g. HD 209458b; 
IFortnev et al.l uOOSl ) . We find that this planet's emis- 
sion spectrum is reasonably well-described by a simple 
1840 K blackbody function, and further note that it de- 
viates from a blackbody in ways that are more character- 
istic of an inverted atmosphere, with a fiux that is lower 
than predicted in the 3.6 /im band and higher in the 
4.5 /im band (see 21) . For the purposes of this paper we 
therefore place this planet in the inverted class, and pre- 
dict that more generalized models for temperature inver- 
sions (e.g. iBurrows et al.|[2008t iMadhusudhan fc Seageii 
[2009) should provide an improved fit to the observed fea- 
tures, albeit with marginal statistical significance given 
the relatively large uncertainties in the planet's observed 
emission spectrum. 

Unlike the ambiguous case of WASP-2b, CoRoT-2b has 
a well-measured emission spectrum that do es not match 
the predictions of either class of models (jGillon et al.l 
I20ld iDeming et all 120101 ). The 8.0 fim. flux for this 
planet is anomalously low relative to its 4. 5 iim flux (for 
a dis cussion of possible explanations see IDeming et al.l 
l2010f ). and we therefore denote this planet in Fig. |3]and 
[5] with a gray square rather than placing it in either of 
the two standard classes. 



We exclude planets with obser vations in only 
one bandpass (e.g. HD 149026b; [Harrington et al.l 
I2007t iKnutson et all l2009a^ , as well as systems with 
only ground-based obse rvatio ns (e.g. WASP-19b; 
iGibson et al.ll2010l lAnderson et a l. 2010), as models for 
these planets are poorly constrained. We also exclude 
cooler (< 1000 K), core-dominated planets such as GJ 
436b and HAT-P-llb, as these planets may not fit into 
the simple classification scheme used to describe hot 
Jupiters atmospheres (indeed, this appears to be the case 
for GJ 436b; Stevenson et al. 2010; Bcaulieu et al. 2010). 
We classify each of the well-characterized hot Jupiters 
as "inverted" or "non-inverted" according to the best- 
consensus interpretation available in the literature for 
that planet, except where otherwise noted. Fig. |4] shows 
the stellar effective temperatures Toff and Ca II H & 
K line strengths for each of these systems. We find that 
the "non-inverted" atmosphere types are consistently as- 
sociated with the most chromospherically active stars, 
whereas "inverted" atmosphere types are associated with 
quiet stars. 

4. AN EMPIRICAL METHOD FOR CLASSIFYING HOT 
JUPITER EMISSION SPECTRA 

The classification scheme described in iJ3]relies on com- 
parisons between the observed hot Jupiter emission spec- 
tra and ID model atmospheres for these planets. In 
this section we develop a model-independent classifica- 
tion scheme that relies on the relative secondary eclipse 
depths in the 3.6 and 4.5 fj,m IRAC bandpasses to dis- 
tinguish between atmosphere types. We select these two 
bands because they are (a) some of the most widely 
available observations for hot Jupiter atmospheres, (b) 
are usually the most precisely-measured eclipse depths 
available, and (c) provi de the strongest constraints o n 
atmosphere models (e.g. IMadhusudhan fc Seageii l2009f ). 
Because the 3.6 /zm data for TrES-1 are currently un- 
published, we completed a preliminary analysis of the 
secondary eclipse in this channel using the methods de- 
scribed in Kmitson et al.l ([2009a) and derive an eclipse 
depth of 0.085% ± 0.013%, which we use for the analysis 
presented here. 

In order to classify each planet, we take the measured 
planet-star flux ratios in the 3.6 and 4.5 /xm bands and 
fit them with a blackbody model for the planet where 
the temperature is allowed to vary fre ely. We use a 
PHOENIX model atmosphere (Hauschil dt et all 119991 ) 
for the star where we interpolate in Toff and log(g) to 
match the star's observed properties and set the metallic- 
ity equal to zero. We calculate the predicted planet-star 
flux ratio in each of the Spitzer bands using the transmis- 
sion functions provided in the Spitzer Observer's Manual 
and then compare the resulting predictions to the ob- 
served slope across the 3.6 — 4.5 ^m bands. Planets with 
strong inversions have lower fluxes at 3.6 /itm and higher 
fluxes at 4.5 ^m compared to our best-fit blackbody 
function, while planets without inversions show stronger 
emission at 3.6 fim and weaker emission at 4.5 /im. In 
the models described in ^ this difference results from 
the fact that the 4.5 /im band contains both water and 
CO bands while the 3.6 /im band is relatively unaffected 
by these molecules. As these molecules switch from ab- 
sorption to emission, it alters the slope of the planet-star 
ffux ratio across these two bands accordingly. We show 



6 



Knutson et al. 



-4.4 - 
-4.6 - 
-4.8 - 



o 



-5.2 



-5.4 



-5.6 



-0.2 



3.6 



-0.1 0.0 0.1 

■ 4.5 Micron Slope vs. Best-Fit Blacl<bocly (%) 



Figure 5. Empirical index for classifying hot Jupiter emission 
spectra vs Ca II H & K activity index log(i?^j^). We derive our 
index by fitting the 3.6 and 4.5 /im secondary eclipse depths as 
measured in the Spitzer IRAC bands with a blackbody function 
for the planet and take the difference between the measured slope 
across the 3.6 and 4.5 /im bands and the value predicted by the 
best-fit blackbody function. Planets with blackbody emission will 
have an index close to zero, whereas planets with negative val- 
ues are brighter than predicted at 3.6 fim and fainter at 4.5 ^m, 
and planets with positive values have the opposite behavior. Plan- 
ets classified in the literature as having atmospheric temperature 
inversions are shown as red circles, while planets without temper- 
ature inversions are shown as blue stars (see [[3] for a discussion of 
WASP-2b). CoRoT-2b (grey square) has an unusual spectrum that 
is not well-described by either model, and we have excluded XO-3 
from this plot as i ts ac tivity level is ambiguously determined from 
its spectrum (see i|2.1l l. The vertical dashed line at -0.05 delineates 
a distinction between atmosphere types; in reality there is likely a 
continuum of behaviors ranging from the default strong absorption 
(i.e., no inversion), to weak absorption (corresponding to a weak 
inversion), to no absorption or even emission from an increasingly 
strong stratosphere. 



the resulting difference in slopes in Fig. [Sj where "non- 
inverted" atmosphere types have negative slopes and "in- 
verted" atmospheres have positive slopes as compared to 
our best-fit blackbody functions. Planets whose emission 
spectra are consistent with a blackbody (i.e. no signifi- 
cant absorption or emission) have values close to zero in 
this index. 

This empirical classification scheme provides a good 
match to the "inverted/non-inverted" classifications pre- 
sented in the literature, with the four planets with the 
most negative values in this index (HD 189733b, TrES- 
1, TrES-3, and WASP-4) corresponding to the "non- 
inverted" atmosphere types described above. Planets 
with significantly negative slopes across these two bands 
are associated with the most active stars, providing a 
model-independent confirmation of the proposed correla- 
tion between increased stell ar activity and hot Jupi ter at- 
mosphere types. Although iWheatlev et alj (|201Cl[ ) place 
WASP-2b in the non-inverted class of atmospheres, we 
find that our index for this planet has a value of 0.025 ± 
0.039, consistent with a blackbody or slightly inverted 
atmosphere. As noted earlier, CoRoT-2b's anomalous 
emission spectrum is poorly fit by both inverted and non- 
inverted atmosphere models, although it is possible to fit 
the 3.6 and 4.5 fim eclipse depths with a 1790 K black- 
body as shown in Fig. [SJ 



5. THE CONNECTION BETWEEN CHROMOSPHERIC 
ACTIVITY AND UV FLUX 

The observed correlation between hot Jupiter emission 
spectra and stellar activity levels would seem to suggest 
that the increased UV fiux experienced by the planet 
destroys the high-altitude absorber responsible for the 
formation of temperature inversions. In order to put 
this correlation on a more quantitative footing, we ex- 
plore the connection between Ca II H & K emission line 
strengths and FUV fiuxes in G and K stars. We focus 
our study on Lyman a, as this line is generally the single 
strongest feature in the UV spectra of late-type stars; the 
sun emits as much flux in Lyman a as in all other wave- 
lengths combined shortward of 1500 A ()Fontenla et al.l 
ll988HLandsman fc Simonlll993[) . In this exercise we com- 
pare stars of similar spectral type but differing activity 
levels, and we therefore elect to work with the simpler 
'5'hk values instead of log(i?^j^) as our activity measure. 
log(i?yj^) uses an empirical function of i? — to remove 
the trend towards increasing S'hk values for later spec- 
tral types (caused by the lower continuum fluxes in the 
region of the Ca II H & K lines in these stars) , but this 
correction is unnecessary when comparing stars of the 
same spectral type. 

It is well-established that increased emission in the 
Ca II H & K line cores in late- type stars is an indicator of 
nonradiative heating processes in the chromosphere con- 
necte d to the presence of enhance d magnetic fields (e.g. , 
LeightOT] [T959[ INoves et al.l 119841: ISchriiver et"all 119891: 
Bailunas et al.lll985() . Increased stellar activity as mea- 



sured using Mg II H & K has also been observed to cor- 
relate with increased Lyman a and x-ray emission (e.g., 
[Wood et al. 2005), but to the best of our knowledge there 
are no quantitative scaling laws directly relating Ca II H 
& K line strengths and the integrated Lyman a flux as a 
function of spectral type for a large sample of stars. We 
address this issue by selecting a sample of 26 stars with 
effective temperatures between 5000 — 6000 K observed 
in Lyman a by the Hubble Space Telescope STIS instru- 
ment and its precursor GHRS as described in lWood et al.l 
()2005f ) . For each of these stars we obtain a measurement 
of the Ca II H & K emission line strengths either from 
Keck HIRES observations wh en available, or from stan- 
dard ca talogues published in iDuncan et al.l ()1991[ ) and 
iHenrv et al. (1996). We also include the sun in our ta- 
ble, with an avera ge Lyman a flux of 5.6 x 10~^ times 
its bolometric flux (I Woods et al.ll2000HWood et aL| [2005l) 
and a S'hk value of 0.133 (|Livingston et al.ll2007D . We 
then bin the stars into four categories according to ef- 
fective temperature (5000 - 5500 K and 5500 - 6000 K 
bins) and activity level, where active stars are defined 
as those with S'hk> 0.3 for the cooler temperature bin 
and S'hk> 0.2 for the higher temperature bin. We then 
calculate the Lyman a fiuxes correspon ding to e ach bin 
using the values given in Table 3 in Wo od et al.l ([2005), 
where the authors have used the Lyman a line shape in- 
formation to correct for the effects of absorption from 
the ISM and any stellar outflows. The UV spectra used 
by Wood et al. are based on a collection of archival HST 
observations and tend to be biased towards more active 
stars; our active bins have 8 cooler stars and 11 hotter 
stars, while the quiet star bins each contain four stars. 
We find that the average Lyman a fiux as a fraction of 
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the total bolometric flux is (1.1 ± 0.2) x 10 ^ for quiet 
stars between 5000 - 5500 K and (4.2 ± 0.2) x 10^^ for 
active stars in that same temperature range, where the 
quiet stars have an average S'hk of 0.171 and the ac- 
tive stars have an average value of 0.545 and we set 
the uncertainty on the Lyman a flux estimates equal to 
the standard deviation of the Lyman a fluxes for the 
stars in that bin. This range in Shk is comparable to 
the observed difference between HD 189733 (5100 K, 
'5'hk= 0.508) and its closest quiet analogue WASP-2 
(5230 K, S'hk= 0.159). By analogy, then, we expect that 
HD 189733 will show an enhancement in Lyman a flux 
of approximately a factor of four relative to WASP-2. 
Because the two planets in these systems orbit at nearly 
identical distances from their host stars (0.031 A.U. in 
both cases), we expect the incident UV flux at the sur- 
face of HD 189733b to be a factor of four higher than at 
the surface of WASP-2b. 

We repeat this same experiment for the 5500 — 6000 K 
temperature bins and find an average Lyman a fiux of 
(5.0 ± 0.6) X 10~^ times the bolometric flux for quiet 
stars and (3.6 ± 1.1) x 10~^ times the bolometric flux 
for active stars, corresponding to a factor of eight en- 
hancement for an increase in average S'hk from 0.157 to 
0.354. TrES-3 has an effective temperature of 5530 K and 
Shk =0.305, while its closest quiet analogue XO-1 has an 
effective temperature of 5750 K and 5hk=0.168. Taking 
into account the planet TrES-3b's slightly smaller orbital 
distance (0.023 A.U. vs. XO-lb's 0.049 A.U.), we esti- 
mate that it experiences a Lyman a flux approximately 
thirty times higher than that experienced by XO-lb. HD 
209458 is the only planet-hosting st ar with a direct mea- 
surement of its Lyman a flux in Wood et al.l (2005), and 
has a total flux of 6.9 x 10~^ times its bolometric flux, 
consistent with the other quiet stars in the 5500 — 6000 K 
bin. 

WASP-4 (5500 K) falls on the boundary between tem- 
perature bins, but it is indistinguishable from TrES-3 
(5530 K) and so we estimate its Lyman a flux using 
the hotter temperature bin. If we assume the Lyman 
a flux scales linearly with S'hk, WASP-4 (Shk= 0.194) 
would have a Lyman a flux of 1.1 x 10^^ times its bolo- 
metric flux. Because the planet WASP-4b orbits at a 
comparable distance to that of TrES-3b, this means that 
WASP-4b would receive approximately one third the flux 
of this planet, ten times higher than the flux experienced 
by XO-lb. Comparing it to the stars in the cooler tem- 
perature bin, WASP-4b receives approximately half as 
much Lyman a flux as HD 189733b and twice as much 
as WASP-2b. 

According to the empirical index plotted in Fig, [5l 
TrES-lb lies closest to the boundary between inverted 
and non-inverted atmospheres, and this is reflected in its 
UV flux. The star TrES-1 (5300 K) has a Shk of 0.244, 
corresponding to a Lyman a flux of 1.1 x 10~^ times 
its bolometric flux. This flux is comparable to that of 
WASP-4, but the planet TrES-lb orbits at a greater dis- 
tance and as a result the incident flux at the surface of 
the planet is much less, only 30% greater than that re- 
ceived by WASP-2b. 

Hot Jupiter atmospheres have a relatively high absorp- 
tion cross section in Lyman a and we expect that much of 
this incident flux will be absorbed in the upper regions of 



the atmosphere, where relative abundances of molecules 
may differ from local thermal equilibrium by many or- 
ders of magnitude. Because photolysis cross sections and 
return reaction rates are species-speciflc and can encom- 
pass multiple pathways, it is difficult to predict a priori 
the effects that a significantly increased Lyman a flux 
would have on the relative abundance of the unknown 
stratospheric absorber. We suggest that this would be 
best addressed by comprehen sive photochemistry mod- 
els such as those des cribed in Za hnle et al.l ()2009l I2010D 
and lLine. Liang, fc Y ung (2010). 

6. DISCUSSION & CONCLUSIONS 

We measured the Ca II H & K line strengths for 
flfty transiting planet host stars and correlated these 
activity indicators with planet atmosphere type. Plan- 
ets with non-inverted atmospheres are found orbiting 
the most chromospherically active stars, while planets 
with inverted atmospheres are found orbiting quiet stars. 
The observational correlation described above is indepen- 
dent of the inverted/non-inverted classiflcation scheme 
for hot Jupiter atmospheres, as we find the same re- 
sult for an empirical classification scheme based on 3.6 
and 4.5 /xm secondary eclipse depths alone. We es- 
timate the statistical significance of our result using 
a rank order test in which we randomly select four 
stars out of our sample of fifteen systems observed with 
both Keck/HIRES and Spitzer after excluding XO-3 and 
CoRoT-2b and ask whether those four stars have higher 
log(i?yj^)values than the remaining eleven stars. We then 
repeat this test for a million different trials and find that 
we obtain the four largest log(i?HK;)values for our non- 
inverted planets by chance only 0.13% of the time, corre- 
sponding to a 3.4cr result where the significance is limited 
by the small size of our sample. 

We find that active G and K stars likely exhibit an 
enhancement of 4 — 7 in their Lyman a emission rela- 
tive to quiet stars with similar effective temperatures, 
and suggest that the increased UV flux from active stars 
destroys the high-altitude absorber responsible for the 
formation of temperature inversions. This theory is con- 
sistent with the sulfur model proposed by I Zahnle et al.l 
()2009() , but could also be applied to a range of other sce- 
narios. Because Ca II H & K and UV emission are only 
indirectly linked, deflnitive conflrmation of this theory 
will require an estimate of the UV fluxes incident at the 
surfaces of these planets, taking into account both the 
spectral types and activity levels of their host stars as 
well as their orbital distances. Orbital eccentricity may 
also prove to be important, as planets such as XO-3 likely 
exist in a pseudo-synchronous spin state in which inci- 
dent radiation is distributed more evenly over the surface 
than in the tidally locked case. If this is true, eccentric 
planets may require higher UV fluxes in order to sup- 
press their dayside temperature inversions as compared 
to their tidally locked counterparts. We will address this 
issue in more detail in a follow-up paper, in which we 
examine the relationship between spectral type, Ca II H 
& K emission, and Lyman a emission and derive indi- 
vidual UV flux estimates for the transiting planets listed 
in Table □ 

Although Spitzer exhausted the last of its cryogen in 
May 2009, the 3.6 and 4.5 fim channels on IRAC are 
still functioning at full sensitivity and the telescope con- 
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tinues to operate as part of an extended warm mission 
using those two channels. As part of this warm mis- 
sion we are currently executing a program to observe sec- 
ondary eclipses of 19 additional transiting planets not ob- 
served during the cryogenic mission (ES program 60021, 
PI Knutson). When combined with the remaining un- 
published secondary eclipse data from Spitzer's cryogenic 
mission and other ongoing warm mission programs (e.g. 
ES 60028, PI Charbonneau; DDT 60003, PI Harring- 
ton), data for 20 — 30 additional transiting planet sys- 
tems should become available by summer 2011. We pro- 
pose that the empirical scheme described in ij4] can be 
used to provide an initial classification for these new sys- 
tems, while the increasing availability of ground-based 
secondary eclipse detections as described above should 
help to provide better constraints for atmosphere mod- 
els in cases where we are limited to the two shorter- 
wavelength Spitzer bands. Because active stars are less 
common in the sample of transiting planets (perhaps due 
to selection biases among the ground-based photomet- 
ric surveys and the increased difficulty of radial velocity 
follow-up work for such systems), we suggest that the 
analysis of secondary eclipse data for planets orbiting 
active stars be given a high priority within the current 
large ongoing Spitzer programs. 
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